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The hole energy spectrum with the ion of an acceptor impurity in the quantum dot has been calculated using 
the spherical approximation of the multiband Luttinger model. The dependence of the hole energy levels on 
the impurity location in the quantum dot has been studied. The effect of the impurity location on the dipole 
momentum and the oscillator strength has been analyzed. The hole interlevel absorption coefficient has been 
calculated. 
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1. Introduction 

For the last decade numerous theoretical and experimental works related to quantum dots (QD's) 
have been carried out jlfi^. The physical properties of the spherical QD, such as a dipole transition, oscil- 
lator strength, and optical absorption coefficient can significantly depend on the presence of impurities 
in QD's. Today, a lot of theoretical and experimental works are available, where the effect of impurities 
on the QD optical parameters has been studied. Most theoretical works analyzed hydrogenic (donor or 
acceptor) impurities 1 6-^ld\ located in the center of the QD or quantum anti-dot 1 11]. For a donor impurity 
within the effective mass approximation, an exact solution of Schrodinger equation has been obtained. 
As concerns the acceptor impurity, the multiband model of the valency band has been used, such as Lut- 
tinger model (l^[l3fl. In the case of a harsh change of the heterosystem dielectric permittivity, the theory 
has been constructed which makes it possible to regard the effect of polarization charges on hydrogenic 
impurity fill, E3- 

In general, the impurity can be anywhere in the QD or even outside. In this case, the Schrodinger 
equation does not have an exact solution. That is why variational methods are used in calculating the 
ground and exited states of an impuri ty. It was found that the impurity shift from the QD center is caused 
by the splitting of degenerated levels (3. EH- The number of splitting levels is equal to the number of 
magnetic number values. The energy spectrum of an impurity electron has been so far studied in the 
cubic QD using a variational method 11611 . General properties of an impurity in the spherical QD in the 
presence of an electric field were analyzed in fl7] using a plane wave basis. The Stark effect was re- 
searched. An ellipsoidal QD with off-central impurity in the parabolic potential well is considered in 
Impurity eigenvalues were defined by expanding the exact wave function over exact functions of the 
system with a central impurity. The binding energy was presented as a function of the impurity position 
and ellipticity constant. 

Despite a large number of QD works, we have not come across the works in which the off-central 
acceptor impurity was analyzed within the Luttinger model for cubic crystals. Given that the valence 
band of most semiconductors of cubic symmetry is degenerate in the center of the BrUlouin zone, the 
Luttinger model should be employed for an adequate theoretical analysis of hole and acceptor states. 



* E-mail: leshkoroman@meta.ua 



©V.I. Boichuk, R.Ya. Leshko, I.V. Bilynskyi, L.M. Turyanska, 2012 



33702-1 



V.I. Boichuk etal. 



Moreover, the problem of the off-central impurity is importam because the impurity can be located on 
the surface of the QD. It is known that in many experiments nanoparticles were obtained formed on the 
surface nanoclusters y|. Such problems are an important step for a further construction of the theory of 
surface impurity defects on the QD. 

As noted before, impurities can change optical properties, in particular the absorption coefficient, 
which is caused by the hole or electron interlevel transitions. The interest to intraband interlevel tran- 
sitions is caused by the possibility to use those transitions in the construction of terahertz radiation de- 
tectors |S], because the transition energy is in the terahertz range. Particularly, this urges the study of 
interlevel transitions in the QD with impurities jl^Ql- In particular, in jisll the dependence of the ab- 
sorption coefficient on the impurity position in the QD was studied. It was shown that the displacement 
of a donor impurity ion from the QD center causes shift of the absorption band into the low-energy range. 
Although in these works interlevel transitions in the QD with the donor impurity are described quite in 
detail, but optical transitions caused by the off-central acceptor impurity have not been researched. 

Therefore, the work is aimed at 

• calculating the energy spectrum of the off-central acceptor impurity within Lattinger multiband 
model; 

• studying the QD optical parameters (dipole momentum, oscillator strength, absorption coefficient) 
with the off-central acceptor impurity; 

• analyzing the obtained results with respective results of the donor impurity. 



2. Eigenvalues and eigenfunctions of the off-central acceptor impurity 
2.1. Formulation of the problem and its solution 

We consider a spherical QD heterostructure with an acceptor impurity which is located on the dis- 
tance D from the center of the QD. The radius of the QD is a. Let the heterosystem be constructed of 
cubic crystals with the large band gap Eg and strong spin-orbit interaction A. Let those crystals be 4- 
fold degenerated in the A; = 0. Taking this into account, the spherical Luttinger Hamiltonian is given by 

H=Mri + yp'-r(p-J)' + n(f), (2.1) 



where ji, T - 1/5(272-1-373) are Luttinger parameters, p is the momentum operator, J = i]x + j]y+ ^Jz 
the spin operator j - 3/2. The potential energy of the system is given by 

Tl{r)^W[f,D) + U{r). (2.2) 

The interaction between the acceptor ion and the hole has been written in the Coulomb form 

W{r,D)^ 1 —, (2.3) 

4:7lEoE\f- D\ 

where £0 is an electric constant, e is a dielectric permittivity of the QD. In this work we consider a het- 
erosystem with a large band mismatch. That is why the confinement potential has been chosen as the 
infinitely high spherically symmetric potential well 

f/(r) = |°' 't'' (2.4) 
\ 00, a. 



The exact solution of the Schrodinger equation with Hamiltonian 12. It does not exist. But if the potential 
d2.3) is neglected, the exact solution can be obtained (3 



V/,M(r,0,¥')= v/2/+l ^ (-l)'-J+^4(r)^^ _\, \Y„n{B,cp)xmj, (2.5) 



/=/-;• mi mj 



mi mj -M 
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where h^f[f+ l), h^Hl+ 1), HM, hm, hmj are eigenvalues of operators F^, L^, Fz, L^, Jz respectively, 
Xrrij are spin functions, Y; „, are spherical harmonics, | ^-(^ j are 3-j symbols. Based on the general 



function 1 12. St . exact solutions in the spherically symmetric field have been obtained in yj, 126142^ 
three types of states 



for 



f,M 
II 



,/+l/2, 



/+1/2 



f.M y-'^>'-f M'P/.M 



[e,q)], /S3/2 (/ = /-3/2, /+1/2), 



,/+3/2 , 



/S3/2 (/ = /-l/2,/ + 3/2), 
1/2 (/= 1, 2), 



(2.6) 



where <1>| ^ [d,(p] are spinors which correspond to the spin j - 312. For convenience, as a unit of length 

we use the effective Bohr radius (a^ - 0.53£i7i A), and as a unit of energy we use Ry* - 13.6/(c^7i) eV 
which represents the effective Rydberg energy. In this system of units, radial wave functions have the 
form: 



R'2 



f 



+ A2 



-C-i 



(2.7) 



where / - I, II (the set of numbers determines the type of states), i?J - ji?-^^^'^, Ry j- 

f f— 3/2 f— 1/21 

^2 ~ 1^/ ' ^/ J ' respectively and coefficients have been written with the use of 6-j symbols 



/+3/2[ 



3/2 3/2 



2/(2Z + I)(2/ + 2) 
(2Z + 3)(2/-I) 



Ci = C2(/,/) = MV^(-l) 



3/2+/+/ 



1 + 2 
3/2 



I 1 
3/2 J 

2j 

M= — ■ 

n 



(/+ 1) [1 + 2] 
21 + 1 



ii is the spherical Bessel function of the first type. The solution for the third type of states can be rep- 
resented by the one spherical Bessel function of the first type. Based on the boundary conditions for an 
infinitely high potential well, the dispersion equation has been derived, from which the parameter Xf,i,n 
has been obtained. Here, n is the number of the solution of dispersion equation for other fixed quantum 
numbers. Therefore, the energy of the hole is written in the form 



=!/,Z,n - 



Xf,l,n 



(2.8) 



Since the energy depends on two more quantum numbers, the function <2.6) can be defined 



„iii 



^/M n /' ^/M ~ ^fM n v ^f^M ~ ^fM-n V ''^^^ °^ presence of an impurity in the QD, the Ritz 

variational method was used to determine the energy of the system. To construct the variational function, 
the wave function of the task without the impurity <2.6t - j277t was used. From boundary conditions for 
the function <2.7K A2 was expressed by A\ 



f,M 



A2^-Ai- 



Jl+2 



Jl+2 



(ly/iw^) 

(iV^) 
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Then, i2J) can be written as: 



-C2 



I 



-A 



ii+2 \x 



C2 



-4i 



7;+2 [I v^i + m)] 

Using )2.9> the ground state variational function (first type of states) was written in the form 



(2.9) 



= >li I Qi Si (r) J' (0, (p) + Wi d (r) O^^^'' (0, ] ( 



-ai 



Vr'^+D'^-2rDcose 



(2.10) 



where Qi, Wi are linear variational parameters, ai is the variational parameter, 9 is the angle between 
the direction on the impurity and the hole. The symmetry of the problem makes it possible to choose 
the coordinate system in the way that the axis z passes through the QD centre and through the ion of 
the impurity. Ai can be determined from the normalization condition. For the ground state / - 3/2, 
M= -3/2,-1/2,1/2,3/2. 

The wave function <2.10t in the Schrodinger equation was substituted with the Hamiltonian <2.H . The 
obtained expression was multiplied by the Hermitian conjugate function <2.10t and the final expression 
was integrated by angle variables. As a result, the functional F = F (Qi, Wj, aj, M] was obtained. The func- 
tional was minimized and the ground state energy that depends on | M| was defined. In the case of the 
central impurity D-0, the energy does not depend on | M| . Also, in the case D-0, Ritz variational method 
can be used to calculate the first excited state (second type of states) with the similar wave function: 



'Pii^^ii QiiSii(r)cD 



f.M 



-aur 



which, due to orthogonality spinors <I>^j, ^ [6 , cp) , is orthogonal to the function <2.10) when D = 0. As con- 
cerns the off-central impurity, the wave function cannot be presented in the form similar to <2.10> . Since 
the cos(0) is in the exponent, the functions are not orthogonal. In this case, additional parameters should 
be included for orthogonalization of the functions. Such actions may complicate the minimization of the 
functional. In this regard, we expressed the off-central impurity wave function by the linear expansion 
over the wave function without impurity: 

^>^Y.^iV°' (2.11) 



where y/ . 
tion 



is the wave function of the hole in the QD [W{f,D) - 0] which satisfies the Schrodinger equa- 



The solution of this equation is of the form 1 12. 10> . The index / denotes the whole set of quantum num- 
bers that characterize the state of the hole. After substituting <2.11) with Hamiltonian I I2.H into the 
Schrodinger equation, the linear system of equations was obtained regarding c;: 



E[{Ei-E)Sjt + Wji]ci^O, 



(2.12) 



where Wji is the matrix element of the potential energy <2.3> on functions of the problem without impu- 
rity. The energy spectrum of the acceptor impurity and coefficients c,- were obtained by solving the linear 
system of equations <2.12) and the normalization condition X \ci\^ = 1. 



2.2. Analysis of the spectrum of the acceptor impurity 

Specific calculations were made for CdSe QD. We take the following basic parameters of the crystal 
that form the QD: Eg = 1.841 eV, A = 0.420 eV (H), ji = 1.66, j = 0.41 |30], e = 9.53 I'sT]. Initially the 
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Figure 1. The energy of the central acceptor impurity in the QD. 

energy of the central acceptor impurity D - was determined using the Ritz variational method. Calcu- 
lation results of the energy as the function of the QD radius are presented in figure[T] 

In figure[l] there are two curves which correspond to the ground state (first type of states / - 3/2) - 
curve 1, and to the first excited state (second type of states f - 3/2) - curve 2. Horizontal dashed lines 
denote the energies of the acceptor impurity in the bulk crystal CdSe 1 13]. Those energies saturate at small 
QD radii (a < 40 A). This is due to a small effective Bohr radius (a^ - 8.38 A). As expected, there remains 
a degeneracy by the quantum number M. 




0.0 0.2 0.4 0.6 0.8 1.0 

D/a 

Figure 2. The energy of the off-central acceptor impurity in the QD as a function of the impurity location. 
The QD radius is a = 25 A. 

The calculations of the energy of the off-central acceptor impurity were made using the Ritz vari- 
ational method (only ground state) and by the linear expansion over the wave function without the 
impurity (many states). For a detailed analysis, further calculations E - E{D) were made for the inter- 
mediate QD radius {a - 25 A). Calculation results of the energy of the off-central acceptor impurity are 
presented in figure [2] Solid curves denote the energy of the acceptor impurity which was calculated us- 
ing the method of linear expansion over the wave function without impurity, dashed curves denote the 
same energies that were calculated by the Ritz variational method. As mentioned above, the variational 
functional depends on the quantum number M: F - F(Qi, Wi,ai,M). Therefore, energy of the off-central 
impurity also formally depends on M. But specific calculations showed that the energy of the off-central 
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impurity depends only on | M| . Thereby, if the ion of the impurity shifts from the QD center, the energy 
level of the ground state splits into two levels: with |M| - 3/2 (curve 1') and with |M| - 1/2 (curve 2')- 
The cause of this dependence on |M| can be explained by the violation of the spherical symmetry of the 
problem, and by the preservation of the cylindrical symmetry. A similar splitting of energy levels was ob- 
tained for the donor impurity |15]. However, due to the fact that the ground state of the donor impurity 
does not degenerate in the single-band model, there is no splitting of the ground state energy of the donor 
impurity in contrast to the acceptor impurity 

Apart from the Ritz variational method, there was used a method of linear expansion over the wave 
function without impurity. The first 74 terms were used in calculations. For small DIa, this method with a 
specified number of terms gives the energy which differs by 10 percent from the Ritz method. For a large 
D/a everything changes to the contrary. Since for a small D/a, the symmetry of the problem is close to 
spherically symmetric problem, the contribution of Wji to degenerated states becomes small and in the 
case of D/a - it is equal to zero. In this regard, for small D/a, the number of terms in the expansion 
should be increased four times which dramaticaUy increases the computation time. Therefore, for small 
displacements of the impurity, it is better to use the proposed Ritz variational method, while for a large 
D/a, the method of decomposition turns out to be better for calculations. As an example, for D/a we got 
the following results: 

^l,_3/2 = 0.6091161/^^/2 _3;2.i - 0.644407V'^';2,-3/2;l;l -0-261659l/.!5/2,-3/2;2.0 

+ 0.163952t//^/2,-3/2;i,i - 0-2875061/^" 2,-3/2;l,2 - 0-15593i//^/2,-3/2;2,i + ■■■' (2.13) 

^1,3/2 = 0.609116l/A^/2,3/2;l,0 + 0-644407l/.^'/2,3/2;l;l " 0-261659l//^/2,3/2;2,0 

+ 0.163952i/^i/2,3/2;i,i +0-287506i/.ii/2,3/2;i,2 -0-15593i/aI/2,3/2;2,i +■■■■ (2.14) 

The following terms are an order of magnitude smaller. Those two states have the same energy £1,13/21 = 
-0.50588 Ry*. Similar energies for other values of D/a are indicated by a curve 1 (figure[2). States *l'i,_i/2 
and ^1,1/2 have the energy £1,13/21 - -0.32289 Ry*. Thus, the curve 1 (figure[2) denotes the energy £1,13/21 
of degenerated states ^1,-3/2, ^1,3/2; the curve 2 (figureO denotes £1,11/21 of ^i,_i/2, ^1,1/2; the curve 3 
(figureO denotes £2,13/21 of ¥2,-3/2. *I'2,3/2; the curve 4 (figureO denotes £2,11/21 of ^2,-1/2. *I'2,i/2- 

In practice, it is very difficult to get a single isolated QD. In most cases, the set of QD are received which 
can be characterized by the average QD radius and by dispersion. If the average QD radius is 25 A and 
the dispersion is 16%, the splitting by the quantum number \M\ vanishes. However, the energy of states 
1 and 2 does not intersect each other. This conclusion is important assuming that the dipole transitions 
between states with the energies £1,13/21 and £1,11/21 are forbidden, between £1,13/21, £2,13/21 are permitted, 
and between £1,11/21, £2,11/21 are permitted too. These optical transitions are calculated in the section of 
the paper that follows. 



3. The interlevel transition of the hole of the acceptor impurity 

We have considered the case of the QD irradiated by a linearly polarized light along z direction. The 
dipole transition matrix element between two states is calculated using the function <2.11) . However, 
coefficients c/ are the wave function in the i//9 representation, and we wrote the dipole transition matrix 
element in the same representation 

d°ij^ez\j^e(^y/^^\z\y/°). (3.1) 

After determining the matrix elements of the matrix d", the dipole transition matrix element of the 
off-central impurity was defined 

dnm — {Cn) d Crn , (3.2) 

where Cn is the vector which consists of coefficients c,- for the n-state of the acceptor impurity. The 
calculation of the dipole momentum shows that transitions are possible when the quantum number 
M does not change. That is why di,-3/2;2,-3/2 * (^1-312 ^ ^2,-3/2), '^i,3/2;2,3/2 ^ (¥1,3/2 ^ ¥2,3/2), 
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'^l,-l/2;2,-l/2 + CPl,-l/2 ^ ^2,-1/2). ^^l,l/2;2,l/2 + CPl,l/2 ^ ^2,1/2); | ^^l,-3/2;2,-3/2 1 - pi,3/2;2,3/2 1 . 
I |2 I |2 

pi -i/2;2,-i/2| - |rfi,i/2;2,i/2| ■ Graphics of the matrix element of the dipole moment of the interlevel 
transitions are shown in figure[3] The curve 1 represents |(ii,-i/2;2,-i/2|^ - \ d\,\n\2,\i2\ , the curve 2 rep- 

I |2 I |2 

resents |rfi,-3/2;2,-3/2| - |'^i,3/2;2,3/2| ■ In addition, the oscillator strength of interlevel transitions was 
obtained (figure|4) using the formula 

Pl,M^2,M' - g2^2y^ [^2,\M'\ ~£'l,|M|J |mi,M;2,M'| ■ (3.3) 

The linear optical absorption coefficient which is caused by the interlevel optical transition was defined 
based on the expression jiilE^l : 



«1,|M|;2,|M| M - 0) 





|'^1,|M|;2,|M| 


^ + \dl-\M\;2,-\M\ 




£0£ 


(£2,1^^1 - £'i,iMi - fKof + (ftr)2 



(3.4) 



where w is the frequency of the light, /io is magnetic constant, HT is the relaxation rate caused by the 
electron-phonon interaction and some other factors of the scattering, |M| - 1/2,3/2. The electron density 
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in the QD a is chosen on the assumption that the QD has one hole (impurity hole). That is why a = 
3/ [4na^). We assume that the QD is under low temperature, and its surface is ideally spherical. That is 
why HT can be estimated as the energy level width which is caused by the scattering on acoustic phonons. 
If the temperature of the system is T ^20 K, then HT - k^T ^ 1.7 meV, where is Boltzman constant. 



25000 



20000 
'g 15000 
10000 
5000 



20 40 60 80 100 

hoj, meV 

Figure 5. The optical absorption coefficient. Dashed curves represent i/2;2,i/2 , dotted curves represent 
«! 3/2;2,3/2 ■ The total absorption coefficient is denoted by solid curves. 

The dependence of the optical absorption coefficient on the impurity ion position is presented in fig- 
ure [5] The curves of light absorption are of the form of Lorenz curve. That is why the squares under 
those curves are proportional to the oscillator strength of interlevel transitions. However -Fi,3/2— 2,3/2 > 
f'l, 1/2— 2,1/2. then ai,3/2;2,3/2 > cc\,\i2;2,\i2 for all frequencies. It is seen from figure[5]that the displacement 
of the impurity from QD center causes the shift of absorption bands into a lower energy region. This shift 
is caused by a decrease of the distance between energy levels. Similar results were obtained for a donor 
impurity . However, the difference lies in the absence of the ground state of the donor impurity for 
degeneration as compared with the acceptor impurity. In addition, it should be noted that the absorption 
bands that correspond to a transition between the states with energies £1,13/21 £2,13/21 for the inter- 
mediate Dl a are shifted into a high energy region as compared with transitions £1,11/21 ^ £2,11/21- The 
difference between those transition energies explain the fact that in figure[5] for Dl a - 0.5;0.7, there is 
seen a "structure" in a = aicj) dependence. For a large or small Dl a, the function a - a((i)) is of a Lorenz 
form. In real heterostructures, there is observed a dispersion by the QD size. Calculations show that if the 
dispersion is more than 15%, the half-width of absorption bands will be larger. That is why the structure 
cannot be observed. 
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4. Summary 

In the proposed paper, there was made a theoretical study of the energy spectrum of the off-central 
acceptor impurity based on the Luttinger model. It is possible to determine: 

• the dependence of the energy on the acceptor impurity position in the QD; 

• a partial removal of the degeneracy of energy levels by a quantum number |M|; 

• the QD optical absorption coefficient with the off-central impurity shows that the displacement of 
the impurity from the QD center causes the shift of absorption bands into the low energy range. It 
is also proved that some obtained results are qualitatively simUar to the results of the off-central 
donor impurity. 
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Hei^eHTpaiibHa aKi^enropHa AOMiujKa y c(|)epi/mHiM KeaHTOBM 

TOHL^i 

B.I. BOMHyK, PM. /leiUKO, I.B. Bi.ni/IHCbKI/IM, /l.M. TypflHCbKa 

Kac|)eflpa TeopernHHOi cfiiswKn, flporo6wL4bKWM flepxaBHwR neflarori^HMM yHiBepcMTeT iM. IsaHa ctpaHKa, 
By/1. CipuMCbKa, 3, 82100, /lporo6MH, /IbsiscbKa o6n., yKpaiHa 

y paMKax c(|)epM4Horo Ha6;inxeHHfl 6araT030HHOi MOfle/ii /larTiHflxepa npoeefleHO oS^Kcnenm eHeprernHHO- 
ro cneKTpy flipKi/i 3a Ha^BHOCTi ioHa aKi^enropHoT flOMiiuKu y ccfiepi/iHHiPi KBaHTOBiw Toni^i. floc/iiflxeHO 3a;iex- 
HiCTb eHepreiuHHnx piBHiB flipKW Bifl posTaujyBaHHfl flOMiiuKW y HaHOCucreMi. HpoaHa/iisoBaHO Bn/iwB no/io- 
xeHHfl flOMiiuKM Ha flunoyibHMM MOMeHT la cmy oci4ii/inTopa MixpisHeBnx nepexofliB. Bi/iBHaneHO Koect)it4ieHT 
nor/ii/iHaHHn CBiT/ia, 3yM0B;ieHiiM MixpiBHeBUMM nepexoAaMn fl,\pm. 

lOiioMOBi cnoBa: Koe(f)'mleHT norAmaHHa, HeLieHrpajibHa aKifenTopHa AOMiiuKa 
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